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Abstract 

In shear flows turbulence first occurs in the form of localized structures (puffs/spots) surrounded 
by laminar fluid. We here investigate such spatially intermittent flows in a pipe experiment showing 
that turbulent puffs have a well defined interaction distance, which sets their minimum spacing as 
well as the maximum observable turbulent fraction. Two methodologies are employed. Starting 
from a laminar flow, puffs are first created by locally injecting a jet of fluid through the pipe wall. 
When the perturbation is applied periodically at low frequencies, as expected, a regular sequence 
of puffs is observed where the puff spacing is given by the ratio of the mean flow speed to the 
perturbation frequency. At large frequencies however puffs are found to interact and annihilate 
each other. Varying the perturbation frequency an interaction distance is determined which sets 
the highest possible turbulence fraction. This enables us to establish an upper bound for the 
friction factor in the transitional regime which provides a well defined link between the Blasius 
and the Hagen-Poiseuille friction laws. In the second set of experiments, the Reynolds number is 
reduced suddenly from fully turbulent to the intermittent regime. The resulting flow reorganizes 
itself to a sequence of constant size puffs which, unlike in Couette and Taylor Couette flow are 
randomly spaced. The minimum distance between the turbulent patches is identical to the puff 
interaction length. The puff interaction length is found to be in agreement with the wavelength of 
regular stripe and spiral patterns in plane Couette and Taylor-Couette flow. 

PACS numbers: 
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INTRODUCTION 



Recently substantial progress has been made in understanding the dynamics underlying 
turbulence in pipe flow close to onset . A large number of unstable solutions to the 
problem have been identified (mostly travelling waves) and it has been suggested that the 
turbulent dynamics evolve around such unstable states. Although similar coherent states 
have been observed in experiments [3|, an important qualitative difference between the exact 
solutions observed to date and structures within the turbulent flow is that while the former 
ones are spatially periodic the latter are localized in the streamwise direction. Specifically 
at low Reynolds numbers (Re~2000) turbulence only occurs in localized patches (4, 5| also 
referred to as equilibrium puffs js, 7], which are separated by laminar fluid. Only when the 
Reynolds number is further increased turbulent patches begin to spatially expand. A better 
understanding of the processes underlying localisation [8| is an essential requirement to allow 
further theoretical progress. This regime is also of practical interest since the pressure and 
flow rate fluctuations in the laminar-turbulent intermitent flows are very large (markedly 
higher than in the case of a fully turbulent flow). In addition edge states js-ll] separating 
laminar from turbulent flow and which are potentially relevant for the transition process 
remain bounded in axial direction ever at much higher Re. This indicates that even beyond 
the intermittent regime localisation is relevant to some aspects of the flow. 

Localized turbulence can be identified in a number of canonical shear flows, but unlike 
in pipe and channel flows, for Taylor- Couette (TCF) and plane Couette flow (PCF) regular 
spatio temporal structures have been reported to consist of tilted laminar-turbulent stripes. 
In the case of Taylor-Couette flow the localized turbulent stripes are traditionally known 
as turbulent spirals 12J. The existence of a distinct preferred wavelength has led to the 



speculation of a wavelength instability preceding turbulence [13hl6|. In pipes, on the other 
hand, many earlier investigations of the intermittent regime p, |6[ were focused on quanti- 
ties like intermittency factors and turbulence fraction which for observation times typically 
realized in experiments and simulations (~ a few hundred D/U), are strongly dependent on 
initial conditions as well as imperfections, such as disturbances created at the pipe entrance. 
In this paper we quantify the spatial interaction in the intermittent regime and we aim to 
determine an upper bound on turbulence fractions and friction factors, which unlike above 
quantities do not depend on the specifics of the experimental setup. The resulting puff 
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FIG. 1: Schematic of the experimental setup 
interaction distance is then compared to puff spacings which naturally occur after a sudden 

Reynolds number reduction from a fully turbulent flow. Finally the resulting puff spacing 

and interaction distances are compared to typical wavelength in intermittent Couette and 

Taylor Couette flows. 



EXPERIMENTAL SETUP 



The experimental set up consists of a long glass pipe of D=10 ± 0.01 mm diameter. The 
pipe is made from 1 meter long sections joined to form the total length of 6 meters. The 
junctions are carefully machined to avoid disturbances in the flow. The entrance consists 
of a 200 mm long cylinder with a diameter of D = 100 mm and contains 2 meshes at the 
upstream end and a smooth convergence to the 10 mm pipe at the downstream end. In the 
absence of external perturbations flow remains laminar up to Re = 8000. Distilled water is 
used as the working fluid and formation of algae is avoided by adding 300 ppm of Sodium 
Hypochlorite. 

The flow velocity is regulated by altering the height of a reservoir tank connected to the 
pipe. A pump is employed at the pipe end to recirculate the water to the reservoir and 
ensures that the pressure head is held constant. A sketch of the experimental setup is shown 
in figure [TJ The mean flow velocity can be determined by weighing the mass of discharged 
fluid over a period of time (typically 60 seconds). In order to reduce temperature fluctuations 
that alter the liquid viscosity the water is circulated through a heat exchanger before entering 
the pipe. The temperature of the water is continuously monitored by a PT 100 thermal probe 
which shows typical variations of 0.1 K over 3-4 hours. Repeated measurements show that 
Reynolds number variations are less than ± 0.5 % over a 3-4 hour period. When the flow 
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FIG. 2: Flow visualisation images of a turbulent puff at Re = 2000 at 50 D (top), 100 D (middle) 
and 150 D (bottom) from the point of perturbation. Flow direction is from left to right. The puff 
reaches its typical size in less than 50 D from the perurbation point and the size is constant except 
for statistical fluctuations. 

changes from laminar to turbulent, the drag increases and consequently the flow rate in a 
gravity driven pipe is reduced. This undesired effect is minimized by circulating the water 
through a big flow resistance (in our case a tube of small diameter) before entering the pipe. 
Since ~ 96 % of the pressure drop occurs over the entrance resistance, the mean flow is 
virtually unaltered by the increase of drag during transition in the pipe. Detailed flowrate 
measurements shows that in the present setup the Reynolds number decreases by less than 
10 when the flow is changed from laminar to intermittent. If the flow is in a steady state (i.e. 
laminar or intermittent at a fixed turbulent fraction), so either before or shortly after such 
a reduction, instantaneous pressure measurments (Validyne DP45 low pressure transducer) 
show that fluctuations in the flowrate and hence in Re are even much smaller than Ai?e± 1. 
The largest variations in Re are casued by long term temperature drifts (typically take place 
on order of hours) where a change by 0.1 K results in a Reynolds number change of about 
ARe = 6. 

Five hundred D downstream of the pipe entrance, a segment of the pipe is illuminated 
by a thin light sheet (parallel to the pipe axis). Neutrally buoyant anisotropic particles 
(Mearlmaid AA) are used for flow visualization purposes. These particles have the form 
of flat platelets which align with the shear 17|, ll8j and hence turbulent flow can be easily 
distinguished from laminar motion. Turbulence can be triggered by a brief injection of 
a jet of water through a 1 mm hole in the pipe wall. As shown in figure [2] at Reynolds 
number Re ~ 2000 such a disturbance results in a localized turbulent puff which develops 
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its typical shape in less than 50 advective time units. Images were taken with a 10 megapixel 
still camera and capture a pipe segment of about 10 D long. Throughout the paper non 
dimensional units are used for time and length scales unless stated otherwise. Space is given 
in pipe diameters (D), velocity in the mean flow velocities (U) and time in the corresponding 
unit (D/U). 



RESULTS 



Periodic perturbation 

The first set of experiments are concerned with determining the distance at which puffs 
first interact. When puffs are located too close to each other the downstream puff decays 



19l | . To find the exact distance where this interaction starts periodic perturbations are 
applied to the flow at a fixed position and the rate of puff triggering as response of the 
injected disturbance is observed at a downstream distance. 

The perturbation consists of a water jet injected by a solenoid valve, placed 100 D down- 
stream of the pipe entrance, which is run at frequencies ranging from f = 0.2 Hz to f = 2.5 
Hz. The valve opening time, frequency and water discharge are independently controlled. 
The duration of each jet is tj et = 0.11 s ~ 2.5 D/U and the amount of water injected is 
0.15 ml per discharge which corresponded to less than one- hundredth of the pipe discharge 
at Re=2000. At the distance of 400 D downstream the perturbation point, the number 
of induced puffs, iVp U ff s , is counted during 300 s ~ 7500 D/U. For each frequency of the 
input disturbance the experiments are repeated five times to extend the observation period 
to t bs ~ 37500 D/U. In the following the temporal spacing of perturbations is converted 
into an axial spacing of pulses by multiplication with the mean velocity U, assuming that 
perturbations are advected with the mean flow speed. This assumption is supported by a 



previous study (20J where we show that isolated puffs propagate only ~ 5 % slower than the 
mean flow at these Reynolds numbers. 

The ratio of the number of observed puffs to the number of perturbation pulses is shown 
in figure [3](a) as a function of the perturbation period. For large axial spacing (small fre- 
quencies) each perturbation results in a single puff. However once the puff spacing is below 
~ 20-25 D, new generated puffs interact and annihilate each other and the ratio of created 
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FIG. 3: (a) Ratio of number of puffs to number of perturbations vs. perturbation period. 
(fc)Distance between puffs vs. perturbation period. 

puffs to perturbations drops below one. Figure [3]^b) shows the distance between the gen- 
erated puffs in response to the perturbation period. From figure Mjo) it is again evident 
that the rate of puff generation is in synchrony with the rate of injected perturbation for 
large axial spacing and that the puff spacing is identical to the spacing preselected by the 
perturbation. Once puffs are created at too short distances (< 20-25-D) the number of puffs 
observed downstream drops resulting in an increase of the distance between observed puffs. 
The minimum of each curve in figure [3]^b) corresponds to the perturbation frequency gen- 
erating the maximum number of puffs which can sustain themselves in a given length of 
the pipe. The corresponding distance will in the following be referred to as the optimum 
spacing as it provides the maximum turbulent fraction. For times between perturbations 
shorter than the minimum one (high frequencies), the generation of puffs decreases quite 
sharply for the lower Reynolds numbers. For instance at Re = 2000, only one puff each 
1000 D is observed at high perturbation frequencies (l/i=7D/U). At these frequencies per- 
turbed regions do not always merge to form a single puff but instead quite frequently they 
annihilate each other. This decrease of generated puffs is smaller for the higher Reynolds 
numbers where the distance between puffs seems to saturate for large frequencies. 

Next we investigate whether the generation of puffs depends on the intricacies of the 
jet injection method i.e. the water injected per perturbation and the valve opening time. 
The puff spacing is found to remain unchanged when the jet discharge is varied from 0.15 
ml/discharge to 0.25 ml/discharge. The impact of the jet duration on puff spacing has also 
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FIG. 4: Impact of valve opening time on puff spacings. Three different valve timings (0.11, 0.22 & 
0.33 sec) were used. In the x-axis t* ert is the time between pulses taking valve opening time into 
account (see text) 

been investigated. Here the valve opening time is varied from t va i ve = 0.11 s to 0.22 s and 
to 0.33 s for Re=2000. From figure H] it is evident that all curves collapse when plotted as 
a function of the time between perturbation pulses i.e. the time between the end of one 
jet and the beginning of the next one (tp ert = tpert-^vaivc)- This highlights that the relevant 
time is not the period of the perturbation but the time for which the flow is unperturbed. 
Since the minimum of all curves is the same, the optimum puff spacing is independent of 
the details of the perturbation. 

In the following the friction factor is measured for three cases: for a flow without per- 
turbation, when the flow is tripped by an obstacle and for the periodic perturbation de- 
scribed earlier. The friction factor quantifies the friction losses and it is defined in a pipe 
as / = Ap/(l/2pU 2 L/D), where Ap is the pressure drop over a length L and p is the 
fluid density. When the flow is not perturbed, the friction factor closely follows the Hagen- 
Poiseuille law / = 64/ Re. For a fully turbulent flow the friction factor lies on the Blasius 
curve / = 0.3164i?e -0 25 . On the other hand, in the transitional regime the friction factor 
strongly depends on initial conditions as well as experimental imperfections. To demonstrate 
this dependence, friction factors resulting from two different obstacles are shown in figure 
[5]^b)(red squares and grey open circles). For the small obstacle transition sets in at ARe ~ 
1500 larger than for the large obstacle. Depending on the perturbation amplitude and the 
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FIG. 5: (a) Friction factor vs. perturbation period, (b) Friction factor vs. Reynolds number using 
different perturbation types as indicated in the legend. 

observation point any transitional friction factor (i.e points between the Hagen-Poiseuille 
and the Blasius curves) can be observed for Reynolds number Re>2000. 

In order to investigate the dependence of the friction factor on initial conditions we 
study the pressure drop for varying frequencies of the jet perturbation. The flow is again 
periodically perturbed as described earlier and 100 D downstream of the perturbation point 
the pressure drop is measured over a length of 770 D (the pipe length was increased by 5 
meter for this purpose) with a Validyne DP45 pressure sensor. The pressure drop over this 
length lies in the optimal working range of the DP45 sensor resulting in an accuracy better 
than one percent for the pressure readings. Starting from small perturbation frequencies, 
i.e. large spacing, the friction factor is found to first increase with frequency and then 
reaches a maximum at the optimum puff spacing (figure |5^a)). Beyond this point the 
perturbed regions start to interact resulting in a lower density / packing fraction of puffs. 
The maximum packing fraction data provides an upper bound for the friction factor for a 
given Reynolds number, shown by the green triangles in figure E^b). These data therefore 
provide a well defined connection between the Hagen-Poiseuille law and Blasius law. While 
the maximum packing fraction is clearly differentiated for Re<2300 the curves become natter 
when increasing the Reynolds number until the point that the curves are almost flat and 
only a very weak interaction can be detected (Re = 2500). The fading of the maximum 
indicates that a fully turbulent flow state is approached. This is in qualitative agreement 



with recently suggested transition points to fully turbulent flow at Re = 2550 to 2600 
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22| 
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although from our data no criticality can be inferred. 



Reduction experiments 

In the following we would like to investigate if in analogy to TCF and PCF "large- 



wavelength instability" scenario [13l-ll5j can also be observed in pipe flow. In the above 
studies it has been suggested that when the Reynolds number is reduced from the fully 
turbulent regime, an instability occurs which gives rise to periodic laminar-turbulent pat- 
terns. If the same holds for pipe flow, its expected that upon reduction of Re from the fully 
turbulent regime, laminar gaps open up at regular intervals giving rise to a periodic laminar 
turbulent pattern. We therefore initially created a uniform turbulent flow at Re=4500 by 
putting an obstacle into the pipe directly behind the entrance and subsequently reduced the 
Reynolds number to the intermittent regime. By pressure measurements we verified that the 
friction factor lies on the Blasius curve (figure [5b) and hence that the flow is initially turbu- 
lent throughout the pipe. Next a bypass valve was opened allowing half the total discharge 
to flow through a parallel pipe (not shown in figured]). By doing so, the Reynolds number 
in the pipe decreases to the desired value. Care was taken to avoid air pockets or elastic 
components in the setup which would cause damped oscillations in the flow rate during the 
reduction. Flow visualization and monitoring of the outflow angle at the pipe exit showed 
that the flow changes smoothly between both Reynolds numbers during approximately four 
seconds (» 100 D/U). 

Images were captured at 10 Hz with a CCD camera (Matrix Vision Blue Fox 121G) placed 
500 D downstream of the entrance. Figure [6](a) shows typical snapshots for a laminar and 
a turbulent flow. The standard deviation in the turbulent image intensity is greater than 
in the laminar flow which allows for a simple identification of the turbulent regions. The 
laminar turbulent pattern in the pipe was reconstructed using the Taylor's hypothesis at 
the Reynolds numbers of the experiments. Generally it was found that after reduction the 
flow rapidly (<100 D/U) adjusts to the new Reynolds number resulting in a fixed pattern of 
equilibrium puffs separated by laminar flow. Consequently, initial 4 seconds (100 time units) 
of each signal were disregarded to allow the flow to establish the equilibrium pattern allowing 
a measurement time of 400 D/U. Once established, this pattern shows little dynamics on 
the time scale of the experiment for Re<2400. This was confirmed by a second set of 
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FIG. 6: (a) Flow visualization at Re = 2000 in turbulent (upper image) and laminar flow (lower 
image). Images display a pipe segment of 5-6 D long. Flow direction is from left to right. One 
reduction experiment is shown in movie 1 available in the online version of the paper, (b) Standard 
deviation of image intensity vs. reconstructed pipe position at Re = 2100. Summits of the peaks 
(black dots) indicate the position of a turbulent spot. 

measurements where the waiting time was increased to 300 D/U resulting in a remaining 
measurement time of 200 D/U (data not shown). Both sets of measurements (200 D/U and 
400 D/U) resulted in the same size distributions (within experimental error). The dynamics 
of puffs only appear to be relevant on much longer time scales, puff splitting in pipes of 
similar length has been observed (231 for Re > 2450 and the probability of an isolated puff 



24j-|26|. A sample signal 



to decay during 500 time units at Re = 2000 is around 5 ■ 10~ 4 
at Re 2100 is presented in figure Mb). The data has been postprocessed by smoothing the 
original signal. The peaks in figure Mjo) represent the presence of turbulent puffs and the 
intermediate valleys correspond to the laminar regions. The summits of the peaks were 
traced and represent the puff positions. After 600 time units all the fluid which initially had 
been in the fully turbulent state had been advected past the pipe exit and the entire flow is 
laminar. 

For each Re 30 to 50 experiments were performed and this resulted in a sample size of 
about 200 puff spacings. In general, the puff spacing decreased with increasing Reynolds 
number. Typically 4-5 puffs appeared for Re = 2000 whereas the number of puffs for 
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Re = 2400 was 8-10 for the same time interval. From the data we determined the cumulative 
probability P(d) that two adjacent puffs are separated (peak to peak distance) by a distance 
d or larger. For a spatially periodic pattern such cumulative distribution should show a sharp 
drop at the wavelength of the pattern. Instead, as shown in figure [7(a), the data falls on 
exponential tails. The fact that the spacing between the puffs is highly irregular is also seen 
in the example shown in figure [6(b). The exponential distributions ( figure [7(a)) do not 
cross zero but reach P(d)=l at some value d=d m i n > where d m i n is the minimum puff 
spacing. The distribution takes the form of : 



P(d) oc e -s for d > d (1) 

where do is related to the distance from which the exponential distribution applies and 
5 is the characteristic scale of the exponential distribution. Note that 5+do is the mean 
distance between the puffs. The statistical analysis of an exponential distribution which 



only applies from a distance d Q has been performed using techniques described in [26j|. This 
probability function implies that the spacing between two puffs is always larger than do but 
otherwise arbitrary and that it is not influenced by the spacing between puffs further up 
and downstream. Therefore the existence of a preferred spacing or wavelength is ruled out. 

For the lower Reynolds numbers (Re = 1900-2200) the exponential distribution applies 
from the minimum puff distance (d = d m i n ) showing that the appearance of a puff is not 
influenced by any upstream or downstream puff. Note that for Re = 2300 and 2400 there 
is a slight initial rounding indicating that size distributions only become random for puff 
spacings slightly larger than d m i n . The deviations from the exponential distributions at long 
distances on the other hand can be attributed to the finite size of the pipe. 

To rule out any effects of initial Reynolds number on the reduction experiments, a second 
set of experiments was carried out where the Reynolds number was reduced from a turbulent 
flow at Re=6000. Excellent agreement is found with the reduction experiments from Re 
=4500 (see the inset in figure 17(a)). Note that the deviations at the end of the tails carry 
little statistical weight and can be attributed to finite sample size. This indicates that the 
size distribution of the observed flow pattern does not depend on the initial intensity of the 
turbulence but only on the Reynolds number established after reduction. 

The characteristic distance S, from equation [fl clearly decrease with Re (figure [7(b)). As 
expected turbulent spots are more closely spaced with the increase of the Reynolds number. 

11 



(b) 



I 



i 



2000 2200 2400 2600 

Re 

FIG. 7: (a) Probability that puffs are separated by a distance greater than d after reducing the 
Reynolds number from Re = 4500. At the initial Re the flow is fully turbulent. The inset compares 
experiments with two initial Reynolds numbers (black triangles for Re = 4500 and green circles for 
Re = 6000). Note that deviations for d >100 were due to the finite sample size, (b) Characteristic 
distance 8 from equation [1] as function of the Reynolds number. Error bars correspond to 95 % 
confidence interval. 

While the quality of the data does not allow to determine the exact functional dependence 
the characteristic distances appear to approach zero for Re >~ 2500, where maximum puff 
packing should be achieved. This point is very close to the transition to uniform turbulence 
(Re c = 2550 — 2600) and to the Reynolds number where the puff interaction diminishes 
(figure Eta)). 

DISCUSSION & CONCLUSIONS 

The experiments presented in this paper establish the existence of a minimum distance 
between turbulent patches in pipe flow. In the first set of experiments an optimum distance 
corresponding to the densest packing fraction of puffs has been determined with the aid of 
a periodic perturbation. Specifically we measured the maximum achievable friction factor 
in the transition regime as a function of Re. Typically the friction factor in this regime 
strongly depends on initial conditions and varies from one experiment to another. Our data 
gives an upper bound for the sustainable friction factor and provides a well denned link 
between the Hagen-Poiseuille and the Blasius law. In the second set of experiments, where 
no periodicity was imposed externally, a minimum spacing has been observed. In Figure 
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FIG. 8: Comparison of minimum puff spacing (Reduction experiments) and optimum puff spacing 
(Periodic perturbation experiments) with laminar turbulent wavelengths in transition regimes of 
Taylor-Couette and Plane-Couette flow from [2]. The upper x axis represents the Reynolds number 
used for Taylor Couette flow and Plane Couette flow and the lower x axis represents the Reynolds 
number used in our experiments. 

[5] we compare the minimum distance from the reduction experiments with the optimum 
distance from the periodic perturbation experiments as function of the Reynolds number. 
Both data sets overlap indicating that the minimum spacing is a property intrinsic to the 
flow, independent of initial conditions or flow perturbations. The existence of a minimum 



19] where it 



spacing in pipe flow can be explained following the arguments presented in 
has been pointed out that the instability mechanism sustaining a puff relies on an inflection 
point at the rear of the puff. When an upstream puff appears the inflection point behind the 
original puff is distorted. Eventually the downstream puff decays if the distortion is strong 
enough. The nature of the instability at the rear end of the puff has also been discussed by 



27] and 11]. Similar to the arguments given in an earlier study by 28| they emphasize the 
relevance of Kelvin Helmholtz vortices created at the rear interface to the puff sustenance. 

It should be pointed out that the active region of a turbulent puff is approximately 10 D 
long [19]. Further downstream the fluctuations have disappeared and the profile gradually 
recovers the parabolic shape until the preceding puff is encountered. This aspect is similar 



to the observations in stripe patterns in PCF [29] . The turbulent eddies are localised in 
the turbulent stripes which approximately have a width of 10 wall spacings. Between two 
stripes the flow relaminarizes but does not quite reach the fully developed laminar profile. A 
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comparison of the minimum puff spacing in pipes to the wavelength measured in PCF and 



TCF laminar turbulent patterns 



13] surprisingly shows that the wavelength closely matches 



the interaction distance found in pipes when the Reynolds number is rescaled so that the 
intermittent region for the different flows overlap. The wavelength have been calculated 
from the spanwise and streamwise spacings measured by [fjj]. Banded laminar-turbulent 
patterns with similar wavelengths have been reported in TCF, PCF, torsional Couette flow 



and plane Poiseiulle flow 



291 ] . Our measurements extend this similarity to the case of pipe 



flow where the typical length scale does not appear as a wavelength but as a minimum or 
optimum distance. We hope that this observation will stimulate further studies to clarify 
similarities and differences between these flows. 

Our results indicate that in pipes the spacing between adjacent turbulent regions is a 
consequence of nearest neighbour interaction process and that no long range instability as 
suggested for PCF and TCF [3, Q is required. It should also be pointed out that as a 
recent study shows [30|] the spreading and splitting of puffs takes place on time scales much 
larger than those available in typical pipe experiments. It is therefore possible that in the 
asymptotic time limit a more regular puff spacing is reached, which may be closer to the 
maximum packing fraction. If so this flow pattern would result from nearest neighbour 
interaction and not manifest a least stable wavelength of an instability. 

Acknowledgments: We thank Dr. Marc Avila and the reviewers for their valuable 
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